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Abstract 

Using a generalization of Cayley's hyperdeterminant as a new measure of tripartite fermionic 
entanglement we obtain the SLOCC classification of three-fermion systems with six single particle 
states. A special subclass of such three-fermion systems is shown to have the same properties as 
the well-known three-qubit ones. Our results can be presented in a unified way using Freudenthal 
triple systems based on cubic Jordan algebras. For systems with an arbitrary number of fermions 
and single particle states we propose to use the Plucker relations as a sufficient and necessary 
condition of separability. 



I. INTRODUCTION 



The quantification of multipartite entanglement is one of the most important problems 
of quantum information theory. Regarding entanglement as a resource proved to be a useful 
idea producing spectacular applications such as teleportation 1 , cryptography 2 and quantum 
computing 3 , and paved the way for further possible fascinating applications. However, this 
idea immediately leads us also to the need of classifying different types of entanglement 
via suitable entanglement measures. These measures are real-valued functions of quantum 
states trying to quantify the amount of entanglement these states contain. For systems of 
distinguishable constituents characterized by either pure or mixed states on the structure 
of such entanglement measures a great variety of results is available 4,5 ' 6 . However, much 
less is known about the structure of multipartite entanglement measures for systems with 
indistinguishable constituents. For bipartite fermionic and bosonic systems a number of 
useful results exists 7 ' 8 ' 9,10,11,12 . For example for two-fermion systems having 2K single particle 
states a decomposition similar to the Schmidt decomposition was introduced 11 . This Slater 
decomposition uses the concept of Slater rank, i.e. the number of Slater determinants 
occurring in the canonical form of any bipartite fermionic system, for the quantification 
of fermionic entanglement. The simplest nontrivial example occurs for a two fermion system 
with four single particle states. Here the states are characterized by six complex numbers 
P12, Pu, P23, P24, P34 that can be arranged into a 4 x 4 antisymmetric matrix P ab , a,b = 
1,2,3,4. It turns out that states of Slater rank one are the ones for which the Pliicker 
relation 13 

P12P34 — P13P24 + P14P23 = (1) 

holds. From multilinear algebra it is well-known that this condition is a sufficient and 
necessary one for writing P ab in the form: P ab = v a w b — w a v b for some four component 
vectors v and w i.e. in this case P = v Aw, it is a Slater determinant. Such fermionic states 
are called separable. When the quantity in Eq. (1) is different from zero we have states of 
Slater rank two, in this case we have precisely two terms in the Slater decomposition and 
the state is entangled. A useful measure of bipartite entanglement in this case is 11 

< 77 = 8IP12P34 - P13P24 + P14P23I < 1 (2) 
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Notice that the measure satisfies 13 if = 64|Det(P)| hence it is invariant under the action of 
the group SL(A, C) of the form 

P aibl - G ai a 'G b ^P a2b2 , G e SL(A, C), (3) 

where summation for repeated indices is understood. Such transformations form a sub- 
group of SLOCC transformations (stochastic local operations and classical communication) 
introduced by Diir et.al. 14 The SLOCC group is just the group of invertible 4x4 complex 
matrices i.e. GL(A, C). Two states are SLOCC equivalent iff there exists a G G GL(A, C) 
transformation converting one state to the other. Since under a SLOCC transformation 
77 I— > |Det(G)|?7, there are only two SLOCC classes corresponding to the cases 77 7^ and 
r\ — 0. Since 77 = characterizes the separable states there is only one nontrivial SLOCC 
class for two fermions with four single particle states. 

As the first nontrivial case of multipartite fermionic entanglement in this paper we address 
the classification of the simplest of three-fermion systems. By virtue of duality of forms 
/\ 3 C 4 ~ /y 1 C 4 three-fermion systems with four single particle states can be mapped to 
single fermion ones hence the states of such systems are not correlated. (Alternatively, we 
can consider the physically relevant interpretation of a particle-hole transformation as a 
manifestation of this duality 7 .) Similarly /\ 3 C 5 ~ /\ 2 C 5 hence a five dimensional three- 
fermion state can be mapped to a two fermion state. Since the rank of the coefficient matrix 
P a b is always even this case can be related to the four dimensional one 7 and the measure 77 
again can be used. Hence as far as multipartite correlations are concerned these cases are 
not interesting. The first nontrivial case is a three-fermion system with six single particle 
states. From the mathematical point of view these states can be represented by elements of 
/\ 3 C 6 the three-fold antisymmetric tensor product of the six dimensional state space C 6 . 

In this paper we classify different entanglement types of three fermions with six single 
particle states under the SLOCC group GL(6, C). In Section II. we introduce a new tri- 
partite entanglement measure T 12 $ quartic in the 20 amplitudes of our fermionic state that 
we later show to be the natural generalization of the well-known 15 three-tangle r 12 3 playing 
a central role in the classification of three-qubit systems 14,16,17 . Then in the form of two 
Theorems we present our main result: two fermions with six single particle states have four 
SLOCC classes, however only two from these classes represent genuine tripartite entangle- 
ment. The representatives of these classes will be given. Then two further quantities of order 
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three and two in the amplitudes are introduced. Taken together with T123 they provide the 
sufficient and necessary set of quantities to determine which class a given state belongs. We 
illustrate the use of these quantities by a very simple example: two states having the same 
single particle reduced density matrices, however from the tripartite perspective they are 
entangled differently. Our classification has a striking similarity to the well-known SLOCC 
classification of three-qubits. In Section III. we elaborate this point and show that this 
similarity is not a coincidence. We introduce a three-qubit-like state which is a tripartite 
fermionic one with only eight nonvanishing amplitudes. It is shown that our new measure 
7^23 reduces in this case to the three-tangle T123 based on Cayley's hyperdeterminant. This 
analogy with three-qubit states enables a different construction of our basic quantities of 
order four, three, and two related to the ranks of the states appearing in the canonical forms. 

We left the proof of our Theorems to Section IV. The reason for this is that these proofs 
are essentially available in the mathematical literature, however in the somewhat exotic field 
of Freudenthal triple systems and cubic Jordan algebras. The key observation is that there 
is an isomorphism between these Freudenthal triples of a special kind and our three-fermion 
systems. Moreover, this isomorphism lifts equivariantly to an isomorphism between the in- 
variance group of such triples and the SL(6, C) subgroup of the SLOCC group. After this 
observation the SLOCC classification follows immediately from the corresponding classifi- 
cation of the canonical forms of the relevant Freudenthal triples. For convenience we also 
presented in Section IV. a brief summary of the relevant concepts of these mathematical 
structures. The reader interested in the details might consult the references given in this 
section. In Section V. we present our conclusions. Here we also would like to propose the 
use of Plucker relations as a sufficient and necessary condition of separability for fermionic 
systems with an arbitrary number of constituents and single particle states. Using some 
recent mathematical results we emphasize the basic importance of the case of two fermions 
with four single particle states and the associated three-term Plucker relation Eq. (1). In 
some sense the problem of separability for any fermionic system is encoded into an equiv- 
alent two fermion system with four single particle states. Hence as a test for separability 
the measure 77 of Eq. (2) is universal. With some further comments on interesting open 
problems we conclude. 
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II. THE SLOCC CLASSIFICATION OF THREE-FERMION SYSTEMS 



Let us consider three fermions with six single particle states. The Hilbert space for one 
fermion is C 6 , hence the total Hilbert space is H = /\ 3 C 6 i.e. the three-fold antisymmetric 
tensor product of three copies of C 6 . Let us introduce the notation e a A e b A e c for the 
normalized Slater determinant formed from the basis vectors e a , e b , e c , a, b, c = 1, . . . 6 i.e. 
we have 

e a Ae b Ae c = -^(e a ®e b ®e c +e c ®e a ®e b +e%e c ®e a -e c ®e%e a -e a ®e c ®e b -e b ®e a ®e c ). (4) 

We represent a three-fermion state \P) with six single particle states by a three- form P e 
A 3 C 3 as 

1 6 

a,b,c=l 

where the coefficient tensor P a b c is totally antisymmetric hence has 20 independent complex 
components. The condition of normalization yields the further constraint 

|^12 3 | 2 H h |^45 6 | 2 = 1- (6) 

Alternatively our fermionic state \P) can be written as 11 

\p) = E ^jifum (7) 

a,6,c=l 

where f a and fl are fermionic creation and annihilation operators satisfying the usual anti- 
commutation relations. In this case we have J2 abc |"uw| 2 = 1/6, hence w a b c <-> P a b c /V6- 
The group of stochastic local operations and classical communication 14 (SLOCC) is acting 

as 

|P) h-> (G®G®G)\P), GeGL(6,C), (8) 

i.e. we are acting with the same 6x6 complex invertible matrix on our three copies of C 6 . 
This means that for the totally antisymmetric tensor P abc we have 

P i— v C 0,2 C b ' 2 C c ' 2 P (Q \ 

1 a\b\c\ u ai ^61 ^ci 1 aibici- \p ) 

We are interested in finding all the SLOCC equivalence classes of three-fermion states. 
Two states are SLOCC equivalent (hence belonging to the same class) iff their amplitudes 
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P a bc satisfy (9) for some G E GL(6, C). It is convenient to work with the subgroup of 
transformations that have unit determinant. These special SLOCC transformations are 
elements of the group SX(6, C) C GL(6, C). We first determine the equivalence classes 
under SL(6, C) and then find easily how our results modify for the full group GL(6, C). 

Our classification scheme is based on a new tripartite entanglement measure invariant 
under the action of G G SL(6, C) defined in Eq. (9). In order to define this measure we 
reorganize the 20 independent complex amplitudes P a bc into two complex numbers a,/3 and 
two complex 3x3 matrices A and B as follows. As a first step we change our labelling 
convention by using the symbols 1, 2, 3 instead of 4, 5, 6 respectively. The meaning of the 
labels 1, 2, 3 is not changed. Hence for example we can alternatively refer to P^q as Pj^ or 
to P125 as -P 12 2- Now we define 
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Mnemonic: the row index of A transforms to the first, the column index transforms to the 
second and third index of P a b c with the overlined pairs are coming from the corresponding 
complements of the column index in cyclic order. E.g. for A12 the column label 2 is replaced 
by the pair 31. For matrix B we have to form the complement of A, i.e. replacing the 
numbers with their overlined versions. 

The new tripartite entanglement measure for fermionic systems with six single particle 
states we wish to propose is 

< T 123 = |T 123 | < 1 (13) 

where 

T 123 = 4 ([Tr(AB) - a{3} 2 - 4Tr(A^) + 4aDet(A) + 4/3Det(S)) , (14) 
where A* and B* correspond to the regular adjoint matrices for A and B i.e. 

AA i = A^A = Det(A)7, BB i = B^B = Det(S)/, (15) 
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with / the 3x3 identity matrix. The overall factor 4 in Eq. (14) is chosen to ensure 7~i 23 < 1 
for normalized states. 

For the classification of three-fermion states we need to define \P) the dual of our three- 
fermion state \P). For this we define a new tensor P a i> c by defining the corresponding 
quantities (a, (3, A, B) via Eqs. (10-12) 

a = -a 2 (3 + aTr(AB)-2Det(B), (3 = a/3 2 - (3Tr(AB) + 2Det(A) (16) 

A = 2BxA s - 2f3B* - [Tr(AB) - a(3}A, B = -2A xB s + 2a A* + [Tr(AB) - a{3]B. (17) 
Here for two 3x3 matrices we have 

Mx N = (M + Nf -M* -N*. (18) 

Notice that the "state" 



P = l E Pa bc e a Ae b Ae c (19) 

a,b,c=l 

is cubic in the original amplitudes and it does not have to be normalized. It will turn out 
that P a b c has the same transformation properties as P a b c described by Eq. (9). Its role will 
be clarified later. 

Now we can state the main result of this paper. 

Theorem 1: A three-fermion state with six single particle states can be entangled in two 
inequivalent ways. The two classes with genuine tripartite entanglement are characterized 
by T 12 3 7^ 0, and T 12 3 = 0, P ^ 0. States with T 12 3 = 0, P = are either separable or 
biseparable. 

It should be clear that Ti 23 is a complex number but P is a collection of 20 complex 
numbers. The condition P = means that these complex numbers are all zero (i.e. the dual 
state is just the zero state). 

Before making a list of the representatives of each class let us consider an example. Let 
us consider the following two states 18 



* = -^=(V2e 1 Ae 3 Ae 5 + e 2 Ae 4 Ae 6 ), $ = -^(e 1 Ae 2 Ae 3 + e 3 Ae 4 Ae^e 1 Ae 5 Ae 6 ). (20) 
It can be easily shown that the single particle reduced density matrices corresponding to 
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1 



(21) 



these two states are the same 

(2 o\ 

1 
2 
1 
2 
\0 1/ 

hence according to the usual tests of entanglement 19 based on a calculation of the eigenvalues 
of the single particle density matrix, \1> and $ have the same amount of entanglement. 
However, for \1> we have 
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(22) 



and for $ the corresponding quantities are 



-P123 - -P345 - -P3T2 - ^> ^156 - A23 - ~^J~> 

hence a calculation of T i2 3 of Eq. (14) using Eqs. (10-12) shows that 

T 123 W = | T 123 ($)=0. 

Moreover, a short calculation reveals that the dual state $ is of the form 

$ = \*sfl(? A e 3 A e 1 = -^v^e 1 A e 3 A e 5 . 
9 9 



(23) 



(24) 



(25) 



Since this state is not identical to zero $ is neither separable nor biseparable. Hence we 
conclude that the states * and $ are representatives of our two different classes with genuine 
tripartite entanglement. 

Now we get back to the SLOCC classification of three-fermion states with six single 
particle states. We have the following 

Theorem 2. Including the classes of biseparable and separable states we have four disjoint 
SLOCC classes. The representatives of these classes can be brought to the following forms 



P= -(e 1 Ae 2 Ae 3 + e 1 A e 2 A e 3 + e 2 A e 3 A e 1 + e 3 A e 1 Ae 2 ), T 123 (P) ^ (26) 
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P = ^(e 1 Ae 2 Ae 3 + e 1 Ae 2 Ae 3 + e 2 Ae 3 Ae 1 ), T 123 (P) = 0, P ^ (27) 
v3 

P=^(e 1 Ae 2 Ae 3 + e 1 Ae 5 A e 3 ), T 123 (P) = 0, P = (28) 
v2 

P = e 1 Ae 2 Ae 3 , T 123 (P) = 0, P = 0. (29) 

Obviously the last two classes correspond to biseparable and separable states. The represen- 
tative of the second class is very similar to the state $ of Eq. (20). For the representative 
of the first class we prefer the four term form, but we will show later that this class can 
alternatively be represented by a two-term expression (as we also expect from our study with 
the state \1/ of Eq. (20)). Notice also the striking similartity with the well-known SLOCC 
classification obtained for three-qubit states 14 . This is not a coincidence as we will show in 
the next section. 

In order to complete our classification we should find a means of deciding whether a 
state having no genuine tripartite entanglement is separable or biseparable. Indeed, the 
conditions Tl 23 (P) = and P = do not specify whether our state is totally separable or 
merely biseparable. 

Let us call a fermionic state |P) separable if the corresponding form P e /\ 3 C 6 is decom- 
posable i.e. if it can be written as P = uj\ A uj 2 A u 3 for some u>j G C 6 , j = 1, 2, 3. As it is 
well-known in multilinear algebra 20,21 the form P is decomposable if and only if the 

4 

IWP) = ^(-l)*- 1 ^^^ = (30) 
i=i 

Pliicker relations hold, for all A = {ai, a 2 } two and B = {b\, 6 2 , 63, 64} four element subsets of 
the set {1, 2, 3, 4, 5, 6} (or alternatively the one {1, 2, 3, 1, 2, 3}). Here hi means that we have 
to delete 6, from the list &i,&2,&3,&4- The Pliicker relations define a set of quadratic forms 
labelled by all possible subsets A and B compatible with the antisymmetry properties of 
the 20 components P a bc- An excercise in combinatorics 22 shows that the number of possible 
quadratic forms is 45. 

As an example let us consider the state 

Q = he 1 A e 4 A e 5 - e 1 A e 4 A e 6 + e 3 A e 4 A e 5 - e 3 A e 4 A e 6 ), (31) 

with P 145 = P345 = — P 146 = — P 346 = 1/2. The only subset combinations to be checked are 
({14}, {3456}), ({34}, {1456}), ({45}. {1346}), and ({46}, {1345}), all of them give the same 
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associated quadratic form, which is zero. Hence this state is separable, as it has to be since 
fi= I( e i + e 3 ) Ae 4 A(e 5 -e 6 ). 

Now we can complete our classification scheme by identifying the biseparable states as 
the ones having Ti2z{P) — 0,P — 0,U(P) 7^ 0, and the separable ones with T\2^{P) = 
0, P = 0, n(P) = 0. Here n(P) = refers to the vanishing of the relevant 45 quadratic 
combinations of the 20 independent amplitudes P a b c . 



III. AN ANALOGY WITH THREE QUBIT SYSTEMS 

In the previous section we have found a striking similarity between our SLOCC classifica- 
tion of three- fermion states and the one for three-qubit states 14 . Now, by studying a special 
subset of three-fermion systems we show that our classification can indeed be regarded as a 
generalization of the well-known results obtained for three qubits. The subsystem we wish to 
study has merely 8 complex amplitudes. In the notation of Eqs. (10-12) these nonvanishing 
amplitudes are arranged as 
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Due to the antisymmetry properties of the tensor P a b c we can arrange all these ampli- 
tudes to have the 1 or 1 in the first, the 2 and 2 in the second, and the 3 and 3 
in third position. Hence we have a state with a collection of 8 complex amplitudes 
(Pi23,Pi23> A23, Pi23>PT23> P T23> P T23> P T23)- Let us denote this new state by \V) and the 
associated 3-form by 



V = P^e 1 A e 2 A e 3 + P 12l e l A e 2 A e 3 H h Pj^e 1 A e 2 A e 3 . 



(33) 



Let us compare this with the usual expression for a three-qubit state 



ip = V'oooe <8> e° <g> e° + ^ooie u <8> e u <g> e 1 H h ^me 1 ® e 1 ® e 



o <> o, A 



,i ^ „i ^> „i 



(34) 



or in the usual notation of quantum information theory 



^ooo|000) + Vooi|001) + • • • + | 111). 



(35) 
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We see that if the indices 1,2,3 refer to subsystem labels of some fictious system and the 
lack of overbar corresponds to and an overbar corresponds to 1 we have a mapping between 
the three-qubit states and our special three fermion states. 

Let us work out how T 123 (P) looks like. In order to make expressions more transparent 
by an abuse of notation we apply the instructive three-qubit labelling, hence we have Pi 23 = 
Pooo, -P123 = -P001 e.t.c. A further simplification can be obtained by reverting to decimal 
notation, i.e. Pi 23 = Po, P 123 = Pi, • • • , P123 = Pi- Using this notation the final expression 
for T 123 (P) is 

T 123 = AD{V) (36) 

where 

D(V) = (P P 7 ) 2 + (PP 6 ) 2 + (P2P5) 2 + (P3PO 2 - 2(P P 7 )[(P 1 P 6 ) + (P 2 P 5 ) + (P3P4)] 
- 2[(P 1 P 6 )(P 2 P 5 ) + (P 2 P 5 )(P 3 P 4 ) + (P^mPe)] 

+ 4P P 3 P 5 P 6 + 4P 7 P 4 P 2 P! (37) 

is Cayley's hyperdeterminant 23 ' 24 . It is related to the three-tangle 15 r 123 the canonical mea- 
sure of tripartite entanglement as ri 23 = 4\D(V)\. Hence for a normalized special three- 
fermion state we have 

< T 123 (P) = r 123 (P) < 1. (38) 

Now it is easy to understand the similarity between our classification as presented by 
Theorem 2. and the usual one for three-qubit systems. Mapping the representative states of 
Eqs. (26-29) of Theorem 2. to a corresponding three qubit one we get the four possibilities 

i(|000) + |011) + |101) + |110)) (39) 
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-L(|000> + |011> + |101» (40) 

^(1000) + |011)) (41) 
|000). (42) 



It is easy to show that 



-(|000) + |011) + |101) + |110)) = (if<g>if <g> if) -^=(|000) + 1 111)) = H®H®H\GHZ), (43) 
2 V 2 
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and 

-U|000) + |011) + |101)) = (/(g)/(8)X)4=(|001) + |010) + |100)) = (I®I®X)\W). (44) 
v 3 v 3 

where if and X are the usual Hadamard and bit flip gates 

V I ' I ■ ( 15) 





Hence, these states are local unitary (hence also SLOCC) equivalent to the usual GHZ and 
W states 14 . Notice also that since 24 

D((Gi ®G 2 ® G 3 )^) = Det(G 1 ) 2 Det(G 2 ) 2 Det(G'3) 2 J D(V'), G u G 2 , G 3 e GL(2, C) (46) 

none of these transformations changes the value of Cayley's hyperdeterminant. 

Now in order to complete our demonstration that the three-qubit SLOCC classification 
is naturally embedded into the one of Theorem 2. we have to also show that the SLOCC 
group GL(2, C)® 3 is indeed embedded into our GL(6, C). It is obvious that this embedding 
is effected by looking at that subgroup of GL(6, C) that leaves invariant the special form of 
the state V of Eq. (33). Such states are clearly the ones leaving the subspaces 11, 22 and 
33 invariant. The reduction of an element of GL(6, C) hence contains three 2x2 blocks of 
GL(2,C) transformations Gi,G 2 and G 3 . By virtue of the antisymmetry property of the 
tensor P abc it is now easy to see that the action of Eq. (9) gives rise to the usual one of the 
form Gi <g) G 2 <S> G 3 . 

This embedding of the three-qubit system into our three-fermion one is also useful to find 
an alternative expression for our new tripartite entanglement measure 7[ 23 . First recall the 
alternative expression 15 for Cayley's hyperdeterminant 

DW = -ie AlA3 e A2A4 e BlB2 e ClC2 e B3B4 e C3C ^A lBl c 1 V'A 2B2 c 2 ^3i J 3C3^ 4 s 4 c 4 (47) 

where A ± , . . . C 3 = 0,1. Proceeding by analogy the relevant expression we have found for 
T 123 is 

T -1 ;J_ Cu6lCia3&2C2,a2&3C3a4&4C4 p p p p I AQ\ 

± 123 c t 1 aibici 1 C12D2C2 1 a^Ozes 1 a^b^c^i \^°) 

where now a±, . . . c 4 = 1, 2, . . . 6. Notice that this expression can be written in the form 

T 123 = - l -t ahcdef P abc P de}) (49) 
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where P is the dual state introduced in Eq.(16)-(17), with its new form 

p ^ klmk'l'm' p p p ( 

t^abc — r a i m f khc F k i Vm i . [00) 

Notice that Eq.(49) can be used to define the antisymmetric (symplectic) form 

{., •} : (P, Q) {P, Q} = ^e abcdef P abc Q def E C. (51) 

It should be clear that if P abc and Q abc are transforming according to Eq. (9) the symplectic 
form is invariant under the SL(6, C) subgroup of the SLOCC group, and transforms by 
picking up a factor proportional to the determinant for the full SLOCC group. 

Now we can neatly summarize the quantities and the role they are playing in our SLOCC 
classification. We need three quantities, of order four, three and two in the amplitudes P a bc- 
They are given by Eqs. (48), (50) and (30). The four different SLOCC canonical forms have 
four, three, two and one terms (see Eqs. (26-29)). We call these classes of rank four, three, 
two and one respectively. The rank equals four iff T i2 3 ^ 0. The rank is less than or equal 
to three iff T i2 3 = 0, less than or equal to two iff P = 0. Finally the Plucker relation gives 
the result that the rank is less than or equal to one iff H(P) = 0. 

IV. CUBIC JORDAN ALGEBRAS AND FREUDENTHAL TRIPLES 

As we have already discussed in the Introduction the proof of both of our theorems is 
available in the mathematics literature. However, these results are scattered in the exotic 
domain of mathematics of Freudenthal triple systems and cubic Jordan algebras, concepts 
that have not made their debut to quantum information theory yet. The only notable 
exception where these algebraic structures play some role is the current research topic called 
"black hole analogy" where mathematical connections between stringy black hole solutions 
and the theory of quantum entanglement have been established 25 ' 26 . Luckily in order to 
understand the basic correspondence between such algebraic constructs and our fermionic 
systems one does not have to dwell deep into the subject. Here we merely streamline the 
basic ideas of the proof, the interested reader should consult the literature 27 . 

A Jordan algebra J over a field F (we have the complex numbers in our mind) is a vector 
space V over F with a bilinear product o [Jordan product) satisfying the axioms 

A o B = B o A, A 2 o(AoB) = Ao(A 2 oB), A,BeJ, (52) 

13 



here A 2 = A o A. The product is commutative by definition but it does not have to satisfy 
the associative law. The only example that we need in this paper is the trivial one of the 
Jordan algebra of 3 x 3 matrices with complex elements denoted by M(3, C). The Jordan 
product in this case is defined as 

AoB = -(AB + BA), A,BeM(3,C). (53) 

Here AB refers to the usual (associative) product of the relevant matrices. We will be 
interested in the so called cubic Jordan algebras, the ones in which every element satisfies a 
cubic polynomial equation. In our case M(3, C) is obviously a cubic Jordan algebra since 
by Cayley-Hamilton we have 

A 3 - Tr(A)A 2 + -(Tr(A) 2 - Tr(A 2 ))A - Det(A)I = 0. (54) 

In M(3, C) regarded as a Jordan algebra Det(A) is called the (cubic) norm of A denoted 
also by N(A). Moreover, a bilinear form (•, •) : M(3, C) x M(3, C) — > C can also be defined 
by 

(A,B) = Ti(AoB) =Tr(AB). (55) 
One can uniquely define the quadratic adjoint/ sharp map by 

{A\B) = 3N(A,A,B), (56) 

using the linearization of the norm 27 . For our case it turns out that A$ is just the usual 
adjoint (transposed cofactor) matrix familiar from Eq. (15), and the map of Eq. (18) is just 
the linearization of the sharp map. 

Note that in the general case we can construct cubic Jordan algebras via the so called 
Springer construction 27 . In this case one starts with a vector space V with a cubic form 
N : V — > F and a special point c called the base point (in our special case it is just the 
identity matrix I). Then via linearization of N one defines suitable linear, quadratic, bilinear 
and trace bilinear maps that give rise to the definition of the sharp map and its linearization. 
If the trace bilinear form is nondegenerate, and the sharp map satisfies (A tt ) tt = N(A)A then 
we have a Jordan cubic. Then it is proved that every Jordan cubic gives rise to a Jordan 
algebra with unit 1 = c. The Jordan product is given by an explicit formula in terms of 
the linearization of the sharp map, the linear and the bilinear maps 27 . In the following we 
do not need this general construction of cubic Jordan algebras however, it is important to 
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bear in mind that the constructions yielding the canonical forms we are going to describe 
are valid even in this general case. 

Next we define the structure group of the cubic Jordan algebra J as the set of invertible 
F linear transformations g of the vector space J which preserve the norm up to a scalar 
A G F which depends on g only, 

Str(,7) = {ge GL(J)\N(gA) = X(g)N(A),Ae J}. (57) 

For M(3, C) the structure group Str(j7") is generated by transformations of the form 

h: A\-+ A1AA.J 1 , Ai, A 2 G GL(3, C), AeM(3,C), (58) 

and t : A\-^ A T where A T refers to the transpose of A. We denote by Str Q (J") the component 
connected to the identity of Str(J") generated by the transformations h of Eq. (58). 

Now we can define a Freudenthal triple system 28,29,30 . This is a vector space M = M(J) 
constructed from the cubic Jordan algebra in the following way 

M(J) = F©F© J® J. (59) 

Obviously dinuM = 2 + 2dim J . In our case we have 

M = C©C©M(3,C) ©M(3,C), (60) 

with complex dimension 2x9 + 2 = 20. An element of Ai can be written as 



x=\ , a,[3eC, A,BeM(3,C). (61) 




Notice that the quantity x can alternatively be used as a shorthand notation for our quan- 
tities introduced in Eqs. (10-12) related to the 20 amplitudes of our fermionic states. 

On M. there are two important extra structures: a skew-symmetric (symplectic) bilinear 
form {•, •} : M. x M. — > F, and a quartic form q : M. — > M defined by 

{x,y} = a6-P>y + (A,D)-(B,C), x=[ aA ), y=\ lC ), (62) 
q(x) = 2 ((A, B) - a(3) 2 - 8(A", B") + 8aN{A) + 8/3N(B). (63) 
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After recalling that for our cubic Jordan algebra M(3, C) we have N(A) = Det(A) and 
(A, B) = Tr(AB) we see that the tripartite entanglement measure of Eqs. (13-14) is related 
to this quartic form as 

T 123 = 2\q(x)\, (64) 

with x given by Eq.(61) and Eqs. (10-12). Moreover, if we associate to the pair (x,y) 
occurring in Eq.(62) the one (P,Q) of three-fermion states we get 

{x,y} = ±e abaUf P abc Q def , (65) 

which is just the symplectic form of Eq. (51). Recalling the definition of the dual three- 
fermion state Eq. (50) we see that Eq. (48) can be written as 



T123 = -2{x,x} = 2{x,x} = 2q(x), x=\_ I, (66) 




where for the definition of the quantities a,/3,A,B see Eqs. (16-17). In the theory of 
Freudenthal triples x (our dual fermionic state) corresponds to the so called trilinear map 28 
T: MxMxM^F related to the quartic form as q(x) = {T(x, x, x),x}. 

The invariance group of the Freudenthal triple Inv(Af) is the group of invertible F linear 
transformations preserving the symplectic and quartic forms, i.e. 

{gx,gy} = {x,y}, q(gx)=q(x), g G lnv(M), i,t/G M. (67) 

The structure of this group has been studied for example by Brown 28 . It was shown that 
liw(M) is generated by elements of three basic types. We give these generators for the 
case interesting to us i.e. M(J) where J = M(3, C). The component connected to the 
identity lnv (M) of lnv(M) is generated by a (A), 7r(A) and g(A 1 ,A 2 ) where A G M(3, C) 
and Ai, A2 G GL(3, C). The action of these transformations on x G M. takes the following 
form 28 ' 31 




a + (B, A) + (A, A») + PN(A) A + f3A 
B + A x A + /3A* 



, AgM(3,C) (68) 



a A + B x A + «A tt 

B + aA (3 + (A, A) + (B, A*) + aN(A) 



(69) 
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e Ai,A 2 : - Dct(Al) DetfA ) > Al ' A2 G 3, C . (70) 

The total group Inv(A4) is obtained by including in £> also the discrete transformation 
x 1— > 2/ by transforming merely A and by taking their transpose. (See the discussion on 
the structure group of J following Eq. (58)). 

The key theorem for the SLOCC classification of our fermionic systems has been proved 
in the nice paper of Krutelevich 31 . It states that every element of M. is Inv(.M) equivalent 
to one of the following "canonical" forms 

1 diag{0, 0, 0}\ / 1 diag{l, 0, 0} \ / 1 diag{l, 1, 0} \ / 1 diag{l, 1, k} 

/ ' lo / ' \o / ' \o 

'(71) 

where k G C, k 7^ 0. The four cases correspond to the ones based on the concept of rank 
we introduced at the end of Section III. Notice also that using the correspondence given 
by Eqs. (10-12) this classification nearly gives our classification of Theorem2. The only 
subtlety arising is that the rank four case gives an infinity of subclasses labelled by the 
nonzero complex number k. These unnormalized states are of the form 

e 1 A e 2 A e 3 + e 1 A e 2 A e 3 + e 1 A e 2 A e 3 + ke 1 A e 2 A e 3 . (72) 

However, we have not used the full SLOCC group yet. The canonical forms obtained 
for M. use the group Inv (.M) which turns out to be isomorphic to SX(6, C) modulo its 
center 31 . Hence we still have the freedom to rescale our states by using the full SLOCC 
group GL(6, C). From Eq. (14) we see that for our state of Eq. (72) T i2 3 = 16/c, moreover 
from the alternative expression of Eq. (49) we see that T123 picks up a factor corresponding 
to the determinant of the transformation. Hence we can use this extra freedom to achieve 
7i23 = 1 and the canonical form of Theorem 2. 

In order to make the correspondence between Freudenthal systems based on the cubic 
Jordan algebra M(3, C) and fermionic systems with six single particle states precise we have 
to also describe the correspondence between their relevant invariance groups, i.e. Inv (.M) 
and the SLOCC subgroup SL(Q, C). Let us define u 6 = e 27 ™ /6 and u 3 = e 27 ™ /3 the sixth and 
third roots of unity. Then SL(6, C) clearly has a center ujqIq, where Iq is the six dimensional 
identity matrix. Moreover, SL(6, C) transformations of the form 

\P) i-> (w 3 I 6 ) ® (w 3 I 6 ) <g> (w 3 I 6 )\P) (73) 
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leave the state \P) invariant. Hence SL(6, C)/c<j 3 / 6 acts on H = /\ 3 C 6 faithfully. Now the 
dictionary of Eqs. (10-12) provides an isomorphism between the 20 complex dimensional 
vector spaces M. and H. Let us denote this isomorphism by / : M. — > H. Now what we need 
is also an associated isomorphism F of groups F : Inv (.M) — > SL(6, C)/c<j 3 / 6 satisfying 

f(gx) = F(g)-f(x), gelnv (M), x G M. (74) 

Since Inv (.M) is generated by three different classes of elements, we have to give the image 
of these generators under F satisfying Eq. (74). One can check that the relevant map is 31 

'l A ,N 



F : a(A) 



/ 
A / 



F : tt(A') 



/ 



A, A' G M(3,C) 



(75) 



F : e(Ai,A 2 ) i-> AiA 2 



Ai/Det(Ai) 







A 2 /Det(A 2 ) 

where A 1; A 2 G GL(6,C). 

Let us now consider the special case when 



A 3 - = Det(Aj), j — 1,2 (76) 



( [i x X 



A 



n 2 

o o /i 3 y 



^! 0^ 



A' 



^ X 



A! = AT 1 



A 2 
A 3 



(77) 



V2 

o o z/ 3 y 

where fx 1 ,...u 3 G C, and Ai,A 2 ,A 3 G C — {0}. The transformations associated to these 
parameters clearly map a three-qubit-like state of type \V) of Eq. (33) to the same type. 
These transformations leave invariant the subspaces (11), (22), and (33). The 2x2 matrices 
operating on these subspaces are of the form 

1 o\ (l z/ 1j2i3 \ / A li2>3 

^1,2,3 ij ' \o i ) ' \ o Arj )3/ 

These generate three copies of the group SL(2, C) i.e. the SLOCC subgroup of determinant 
one transformations on three qubits. Writing out explicitly the action of this subgroup of 
SL(6, C) on states of type \V) we recover the usual law 



(78) 



\V) ^ {S 1 ®S 2 ® S 3 )\V), 5i, S 2 , S 3 G SL(2, C). 



(79) 



This sheds some more light on our explanation of the three-qubit-like structure embedded 
in the three fermion system we have found in Section III. (See the discussion following Eq. 
(46)). 
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V. CONCLUSIONS 



In this paper we investigated the entanglement properties of three-fermion systems with 
six single particle states. For such systems we introduced a new measure Tl23 of tripar- 
tite entanglement (Eq.(13) ) depending on the 20 complex amplitudes characterizing our 
fermionic state. This entanglement measure is of quartic order, and can be regarded as 
a generalization of the well-known three-tangle T123 of three-qubit entanglement based on 
Cayley's hyperdeterminant. We also introduced two further quantities of order three and 
two in the amplitudes (see Eqs (16-19,30)). They are the dual fermionic state P and the 
Pliicker relations II (P) . Using these three quantities in concert we managed to obtain the 
SLOCC classification of our three fermion systems. We have four SLOCC classes. Apart 
from the separable and biseparable ones we have two nontrivial classes with tripartite en- 
tanglement. The canonical forms of these classes are given by Eqs. (26-29). These states are 
the representatives of the corresponding four classes. For the number of terms appearing in 
the canonical form we coined the term rank. States with 7i 2 3 7^ are of rank four, the ones 
with 7^23 = have rank at most three, the ones with P=0 at most two, and at last fermionic 
states with II(P) = are of rank one. This notion of rank obviously generalizes the con- 
cept of Slater rank well-known from the corresponding classification of bipartite fermionic 
systems. 

We have found a striking similarity between our SLOCC classification and the corre- 
sponding one obtained for three-qubit systems. This is not a coincidence. By employing a 
special three-qubit-like fermionic state with 8 amplitudes we managed to demonstrate that 
the three-qubit SLOCC classification is naturally incorporated within the fermionic one. By 
restriction to the state with merely 8 amplitudes T123 reduces to T123. This phenomenon is 
similar to the one found by Gittings and Fischer 32 for systems of two fermions with four 
single particle states. In this case it is easy to see that the fermionic measure 77 of Eq. (2) 
reduces to the two-qubit concurrence C. This analogy enabled an alternative construction 
for our quantities of order four , three and two providing additional insight into their struc- 
ture (Eqs. 30, 48, 50). Finally we highlighted the proof of our theorems via introducing the 
reader to the basics of cubic Jordan algebras and Freudenthal triple systems. For the proof 
we referred to existing results in the mathematical literature. 

This unexpected connection between such algebraic constructs and quantum entangle- 
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ment might prove to be useful to obtain a classification of further special entangled systems. 
As an example here we mention the recently studied tripartite entanglement of seven qubits 
used in connection with the £7(7) symmetric black hole entropy formula regarded as an en- 
tanglement measure 25 ' 33 . This entanglement measure is again just the quartic invariant q(x) 
(see Eq.(63)) for a Freudenthal triple system however, now it is based on the cubic Jordan 
algebra of 3 x 3 Hermitian matrices with elements taken from split octonions. There is a 
truncation of this formula related to cubic Jordan algebras based on the split quaternions, 
with the correponding entanglement regarded as a truncation of this unusual type of tri- 
partite entanglement. Our results on three-fermion systems fit naturally into this scheme. 
Our cubic Jordan algebra M(3, C) can be shown to be isomorphic to the one based on the 
split complex numbers i.e. the binarions 27 . Moreover, we have already seen that three-qubit 
systems can be regarded as a convenient truncation of this case. In summary all these cases 
of special entangled systems fit nicely into the theory of Freudenthal triple systems based 
on cubic Jordan algebras over split division algebras of complex numbers, quaternions and 
octonions. The question is whether the highly special entangled systems arising in the black 
hole context have any relevance to quantum information theory. 

We would like to emphasize that the nice results that we can obtain for three-fermion 
systems with six single particle states rest on the relationship between these systems and 
the very special structure of Freudenthal systems. So for the general case of multipartite 
fermionic entanglement these structures are not useful. Hence the identification of different 
types of genuine multipartite fermionic entanglement via suitable measures remains a basic 
challenge. However, we would like to point out that as far as the problem of separability 
for fermionic systems with arbitrary number of constituents and single particle states is 
concerned Pliicker relations provide a sufficient and necessary condition of separability 20 ' 21 . 
Hence if we have a /c-fermionic state with n single particle states we can define the /c-form 



Then as usual we call P separable iff P = U\ A uj 2 A ■ ■ • A u k for some Uj G C™. The sufficient 
and necessary condition for this to happen is 

fc+i 

IUfl(P) = W^^W, = 0, (81) 




k 



(80) 
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where A = {a±, a 2 , ... , Qfc-i} and B = {bi, 6 2 , • • • , frfc+i} are fc — 1 and k + 1 element subsets 
of the set {1, 2, . . . , n}, and where the number bj has to be omitted . It is known 22 that we 
do not need to consider all the choices of indices A and B. We have to merely consider the 
elements in A and B in increasing order. If A is contained in B the Pliicker relations are 
identically zero. If we have a single element a G A which is not lying in the intersection of A 
and B we can demand that a < b for all b e B not in the intersection. It is then calculated 
that the number of such subsets is 

1 M ! 

K = ~ A+ h l (lm) (m+ l)!(m + 3)\(k -m- 2)!(n — k — m — 2)V { } 

where M = mm{k,n — k}. So k gives the number of relations to be checked. Moreover, it 
was also shown 22 that one can construct a certain finite set of maps mapping the original 
fc-fermion state with n single particle states to a finite number of two-fermion states with 
four single particle ones. It was shown that the separability of the /c-fermion state is in 
some sense equivalent to the separability of the corresponding two-fermion states. Recall 
the simplicity of the Pliicker relation in this case (see Eq. (1)). So the measure 77 based 
on this relation as far as the separability of /c-fermion states is concerned in some sense 
universal. 

There are a lot of further interesting questions to be addressed. For example it is well- 
known that the three-tangle r 12 3 is an entanglement monotone. What about our newly 
introduced quantity 7i 2 3? There is also the question whether we can generalize our tri- 
partite measure via the usual convex roof construction to obtain a corresponding mixed 
state measure. Moreover, in the original three-qubit setup ti 23 plays the role of the residual 
tangle in the Coffman-Kundu-Wootters relations 15 of distributed entanglement. Can these 
relations be generalized in some sense? We can hopefully address these interesting questions 
in future works. 

VI. ACKNOWLEDGMENT 

One of us (P. L.) would like to thank Professor Werner Scheid for the warm hospitality at 
the Department of Theoretical Physics at the Justus Liebig University of Giessen where a 
part of this work was completed. Financial support from the Orszagos Tudomanyos Kutatasi 



21 



Map (OTKA) (Grants No. T047035, T047041, and T038191) is gratefully acknowledged. 



1 C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres and W. K. Wootters, Phys. Rev. 
Lett. 70, 1895 (1993). 

2 C. H. Bennett and D. P. DiVincenzo, Nature (London) 404, 247 (2000). 

3 M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge 
university Press, Cambridge, england, 2000. 

4 M. B. Plenio and S. Virmani, Quant. Inf. Comp. 7, 1 (2007). 

5 R. Horodevki, P. Horodecki, M. Horodecki and K. Horodecki, arXiv:quant-ph/0702225. 

6 I. Bengtsson, K. Zyczkowsky, Geometry of Quantum states, Cambridge University Press, 2006. 

7 K. Eckert, J. Schliemann, D. Bruss, M. Lewenstein, Ann. Phys. N. Y. 299 88 (2002). 

8 Y.S. Li, B. Zeng, X. S. Liu, and G. L. Long, Phys. Rev. A64, 054302 (2001). 

9 R. Pakauskas and L. You, Phys. rev. A64, 042310 (2001). 

10 G. C. Ghirardi and L. Marinatto, Phys. rev. A70, 012109 (2004), G. C. Ghirardi, L. Marinatto 
and T. Weber, J. Stat. Phys. 108 49 (2002). 

11 J. Schliemann, J. I. Cirac, M. Kus, M. Lewenstein, and D. Loss, Phys. Rev. A64, 022303 (2001), 
J. Schliemann, D. Loss and A. H. MacDonald, Phys. Rev. B63 085311 (2001). 

12 X. Wang and B. C. Sanders, J. Phys. A38 L67 (2005). 

13 P. Levay, Sz. Nagy and J. Pipek, Phys. Rev. A72022302 (2005). 

14 W. Diir, G. Vidal, and J. I. Cirac, Phys. Rev. A62 062314 (2000). 

15 V. Coffman, J. Kundu, and W. K. Wootters, Phys. Rev. A61 052306 (2000). 

16 A. Miyake, Phys. Rev. A67 012108 (2003). 

17 P. Levay Phys. Rev. A71 012334 (2005). 

18 A. Gottlieb and N. J. Mauser, Phys. Rev. Lett. 95 123003 (2005). 

19 See for example the relevant long list of references in Ref.18. 

20 P. Griffiths and J. Harris, Principles of algebraic geometry, Wiley, New York, 1994. 

21 W. V. D. Hodge and D. Pedoe, Methods of algebraic geometry Vol. II., Cambridge University 
Press, Cambridge, 1994. 

22 A. Kasman, K. Pedings, A. Reiszl, and T. Shiota, arXiv:math/0510093. 

23 A. Cayley, Camb. Math. J. 4 193 (1845). 

22 



I. M. Gelfand, M. M. Kapranov, A. V. Zelevinsky, Discriminants, resultants and multidimen- 
sional determinants, Birkhauser Boston 1994. 
P. Levay, Phys. Rev. D75 024024 (2007). 
M. Duff and S. Ferrara, Phys. Rev. D76 124023 (2007). 

K. McCrimmon, A taste of Jordan algebras Universitext. Springer- Verlag New York, 2004 
R. Brown, J. Reine Angew. Math. 236 79 (1969). 
J. Faulkner, Trans. Amer. Math. Soc. 167 49 (1972). 

H. Freudenthal, Indigationes Math. 16 218 (1954), 16 363 (1954), 17 151 (1955), 17 277 (1955). 
S. Krutelevich, arXiv:math/0411104. 

J. R. Gittings and A. J. Fischer, Phys. Rev. A66, 032305 (2002). 
M. Duff and S. Ferrara, Phys. Rev. D76, 025018 (2007). 



23 



